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1.  Introduction 

In  the  theory  of  radiative  transfer  there  are  several  problems 

which  can  be  solved  by  finding  the  solution  H(u)  of  the  non-linear 

integral  equation 

1 
n  1  \  1   _  T  ^  r  g(u)H(u)du 

0 

which  is  called  Chandrasekhar ' s  equation.   The  books  by 
Chandrasekhar  [1]  and  Kourganoff  [2]  contain  discussions  of  this 
important  equation;  and  these  books  also  present  the  contributions 
of  various  mathematicians  who  have  shown  that  (1.1)  can  be  solved 
explicitly  by  using  function  theory  techniques  based  on  analytic 
continuation.   Recently,  (1961 )  C.  Fox  [5]  has  shown  that  (1.1)  can 
be  converted  into  the  linear  equation 

1 


(1.2)  H(x)G(x)  =  1  +  x/  ^i^ 


)du 

X 

0 


where  G(x)  is  known  and  g(x)  is  prescribed.  The  equation  (1.2)  is  a 
singular  Integral  equation  with  a  Cauchy  kernel  and  it  can  be 
solved  for  H(u)  by  using  an  extension  of  Carleman's  method  as  shown 
for  example  in  Muskhelishvili ' s  book  [4]. 

Chandrasekhar ' s  equation  can  be  linearized  by  first  writing 
it  in  the  form 


do) 


r     k (u)du 
xg(x)  =  A<!)^(x)  +(|)^(x)  j    -"^^^ 


and  then  integrating  each  side  of  this  equation  after  it  has  been 
multiplied  by  the  factor  l/{x+i).     This  gives 


ilA) 


r  xg(x) 

J  X  + 


)du 

r 


p    (})^(u)du 
J         u  -  ii 


f 


^   <t)3_(x)dx 


"x+T" 


p"-  (j)^(u)   P-  .t>T(x)dx 


-j   -hi  J   -ir— d- 

0         0 
and  after  multiplying  (lA)   by  (i)-j^(l)  and  using  (1.5)  we  find 


(1.5) 


>;L(e)G^(e)  =  '^  +f 


X-  ii 


This  is  essentially  the  procedure  that  was  used  by  Pox  to  pass  from 
(1.1)  to  (1.2).   It  suggests  the  possibility  of  solving 


(1.6) 


n     (|),  (u)dU 
A^^(x)+^l(x)  j   -inrr^^  h 


(x) 


which  is  both  singular  and  non-linear.   Equation  (1.6),  in  turn^ 
suggests  an  investigation  of  the  more  general  equation 


A^(c)+Mo/i^^=  f(a 


where  C  is  an  interior  point  of  the  simple  smooth  arc  L  which 
connects  the  points  t  and  t,  in  the  complex  r-plane. 

One  of  the  purposes  of  this  paper  is  to  show  in  Section  2 
that  the  equation  I  can  be  solved  explicitly  by  using  elementary 
function  theory  techniques.   It  turns  out  that  the  solution  of  I 
is  in  some  ways  simpler  than  the  solution  of  (1.3).   We  will  also 
be  concerned  with  the  solution  of  some  other  non-linear  equations. 
In  Section  J>   we  show  how  to  solve 


II 


7r2i2(a 


r  c|)(T)dT 

J  T-C 

■-  L 


f(a 


and  Section  ^  is  devoted  to  the   solution  of 

2 


III 


TT^^^iO     + 


/^ 


)d' 


f(a  . 


In  Section  5  we  show  that  there  is  a  connection  between  equations 
I,  II,  III  and  certain  problems  in  potential  theory. 

The  equations  I,  II,  III  may  be  of  interest  for  at  least  two 
reasons.   In  the  first  place,  they  are  of  interest  in  themselves  as 
Cauchy  singular,  non-linear  integral  equations  which  can  be  solved 
explicitly.   In  the  second  place,  they  present  a  formulation  of 
certain  non-linear  boundary  value  problems.   Equation  III,  for 
example,  is  intimately  related  to  a  problem  in  two-dimensional 
potential  theory  which  has  a  number  of  physical  applications.  This 
Is  the  problem  of  finding  a  potential  function  in  a  domain  D  when 
its  normal  derivative  is  prescribed  on  one  part  of  the  boundary  C; 
and  the  magnitude  of  its  gradient  is  given  on  the  remaining  part 
of  C.   In  Section  5  we  show  how  an  explicit  formula  for  the  solution 
of  this  problem  can  be  found. 

The  final  Section  6  is  concerned  with  a  brief  discussion  of 
the  non-linear  system 


L 


and  some  other  systems  which  can  be  linearized  by  the  method 
developed  in  Section  2. 

We  state  here  the  main  conditions  and  assumptions  upon  which 
our  analysis  is  based.   If  t  =  T(t)  is  the  equations  of  the  simple 
smooth  arc  L  directed  from  t  to  t,  ^   t  let  L[t  ,t, ]  denote  the 
set  of  points  T=T(t),  t<tj^t,;  and  let  L(t  ,  t^)  denote  the 
set  L[t^,t,]  minus  the  endpolnts  t  and  t,.   We  assume  for 
simplicity  that  the  unknown  function  (i)(T)  satisfies  a  uniform 
Hftlder  condition  for  any  pair  of  points  in  L(t  ,t, )  and  thus 
guarantee  the  existence  of  the  Cauchy  principal  values  of  the 
integrals  which  appear  in  I,  II,  III  where  C  is  in  L(t  ,t,).  We 
also  assume  that  if  (|)(t)  has  a  singularity  at  an  endpoint  a  of 
L[Tq,t^]  it  is  such  that  £^_^^  (T-a)c|)(T)  =  0.   We  will  say  that  a 
function  with  these  properties  belongs  to  the  class  5V^ .   If  (})(t) 
belongs  to  the  class  9/ ,  it  can  be  shown  that 


/ 


(|)(T)dT 

L 


satisfies  a  uniform  Holder  condition  for  any  pair  of  points  ^,  and 
Cg  in  L(t^,t-);  and 

F(z)  -PA^ 
L 

as  a  function  of  the  complex  variable  z  is  such  that 

The  properties  of  the  prescribed  function  must  match  the 
properties  implied  by  the  representations  on  the  left-hand  side  of 
I,  II,  III.   That  is,  in  accordance  with  our  assumptions  about 


(J)(t),  f(C)  must  satisfy  a  uniform  Holder  condition  on  L(t  ,t-);  and 
although  it  may  have  a  stronger  singularity  than  ^{x)   at  a,   f(0 
must  be  such  that  £^_^^(^-a)  f(0  =  0.  We  will  say  that  under  these 
conditions  f(C)  belongs  to  the  class  ^p. 

The  transformation  2t  =  (t^i-t^q  )v  +  (t^+t^  )  maps  L[t  ,t^] 
into  L[-l,l]  a  simple  smooth  arc  directed  from  t  =  -1  to  t  =  1. 
If  we  also  use  2C,  =    (t,-t  )a>+  (t  +t,  )  we  have 


r         Ux)p  ^  r 

Ut^,x^]  L[-1,1] 


(l>Q(v)dv 


which  shows  that  the  transformation  does  not  change  the  form  of 
equations  I,  II,  III.   Thus  there  is  no  loss  of  generality  if  we 
assume,  as  we  will  hereafter,  that  L  in  I,  II,  III  is  L[-l,l]. 

2.  Equation  I 

We  proceed  to  show  how 

(2.1)  4(o  +  <i)(a/i^^  =  f(a 

L 

can  be  solved  by  using  the  Hardy-Poincare-Bertrand  formula.   This 
formula  states  that  if  ^(t,0  is  suitably  restricted  then 

L       L  L  L 

provided  oo  is  in  L(-l,l).   This  says  that  an  interchange  of  the 
order  of  integration  on  the  left  side  of  (2.2)  leaves  the  residue 
-7r^^(co,a)).   If  <|)(t)  belongs  to  the  class  ':XJ  j^,    defined  in  the 


&fffi   u 


#&*r:i 


•t7-»«  ciaom  (,'■ 


introduction,  the  formula  (2.2)  holds  with  ^(t,C)  =  (1-^  )(|)(0<I>(t). 
The  application  of  (2.2)  to 


1    r     (1-C^)c|)(a  P    <|)(T)dTdC 


l-Oi 


shows,  after  a  little  manipulation,  that 

)dTdC 


r 


-TT^^^iOi)     + 


_/^]^i^[/*-i  • 


Now  if  a  solution  of  (2.1)  exists  when  f(0  belongs  to  '^g^ 
then  (2.1)  and  (2.3)  imply 


(2.4) 


2_   r  (i-r)f(Odc  ^  _iA_   r  (i-r)Mc)dc  .  ^2.2 


^  i^ 


1-co   % 


TTH- 


ir  (b  (o)) 


/ 


^{x)dT 


.l2 


1 
l-(o' 


t2 


T)dT 


(2.5) 


where 


C  -co  M^  \   /       J  t  -  CD 

L 


r  (t)(T)d' 

J  T  -co 


2   [k^ -2A(k^  +  cok^)] 
" 1-co^ 


k,  =/i( 


T)dT 


k^  =  /  T(|>(T)dT  . 

L 


'IT 


AMl'^-' 


The  addition  of 


(2.6)        ±   2Trl   f(co] 


to   (2.5)   gives 


2A7rl   (i)(co)  +27ri   <l)(co)  ^^    iilM] 


(2.7)        -^      r     ^^-^'/^^^^^^   ±    2.i   f(c.) 


±    TTl 


Let  us  introduce  the  function 


[kg-2A(k^  +  cok^)] 

A  +  « 


1-0) 


F^(z)  =  d-z") 


/^ 


)dT 


r      (l-T^)f(T)dT  .g^(  )_^2^(^.^2)^2 

\_y       X  ■■  z  X    o    o 


This  function  is  analytic  for  z  not  on  L[-l,l]  and  it  vanishes  as 
z  — >  CO  .  As  z  approaches  the  point  cd  from  the  positive  and  negative 
sides  of  L,    the  limit  values  of  Fq(z)  are 


F-(a))  =  (1-0)'') 


±  TTl 


*(-'-/  ^^-^' 


2x^/ 


27ri(l-0)2)f(0))  +  2jr  ^^'''^^f';J 


[2A(k^  +^^0^  '^o  +(1-'"^)^^1 


(T)dT 


■;(^    :.-:j\-yl^iii 


but,  as  we  can  see  from  (2.7),  these  limit  values  vanish.  Further- 
more, it  follows  from  the  conditions  vje  have  imposed  on  ^   and  f, 
that  the  behavior  of  F  (z)  in  the  neighborhood  of  an  endpoint  of  L 
must  be  such  that 

(2.8)  £^__^^(l-z)F^(z)  =  0 

Let  n  be  the  boundary  of  a  domain  which  contains  L[-l,l].   If  z  is 
in  the  exterior  of  P  then 


^o(^)-l^i  ^ 


Fo(adc 


Since  (2.8)  holds,  the  path  P  can  be  collapsed  into  the  path  C, 
composed  of  the  upper  and  lower  banks  of  L,  without  changing  the 
value  of  the  Cauchy  integral.   At  these  banks,  however,  the  limit 
values  F"(od)  vanish.  VJe  therefore  conclude  that  Pq(z)  vanishes 
everywhere  in  the  exterior  of  L[-l,l]  and  hence,  if  a  solution  of 
(2.1)  exists,  the  constants  k^  and  k^  must  be  such  that  the  two- 
valued  function 

1 
(2.9)   S(z)  = 


K" 


:l!imMT^  2^(   ^    )  _  ^2  ^  (i_,2),2 


^7^1  "J    ?^ri         -.-n,^^-..^^;  -r.^ 


2 


(A)      =|iL5J|i. 


is  analytic  for  z  not  on  L[-l,l].   We  will  refer  to  this  condition 
as  condition  (A). 


)    ^(':»i) 


■ruaz^^   i 


If  for  a  prescribed  f(0  constants  k^  and  k-j_  exist  such  that 
S(z)  is  analytic  for  z  not  on  L[-l,l],  then,  as  we  can  see  by  taking 
the  limit  values  of  (2.9), 


(2.10) 


(2.11) 


f  ilil|i  +  ^i  i(a+A  =  s+(a 

L 


and,  as  subtraction  of  (2.11)  from  (2.10)  shows,  the  solution  of 
(2.1)  is 


'(a  =  4t  [s+(a  -s-(c)] 


(2.12) 


1 

2Tri 


i-r    ^r     """^ 


)  +  A^ 


1 


,  2A(k,^U,)-k^ 

1  -  c^ 
2     r  (i-T^)f(T)dT  .  2.if(a  +  A^ 

2A(k^  +  Ck^)-k2 


1-C 


For  the  case  in  which  L[-l,l]  coincides  with  the  real  axis  we 
find  that  if  f(|)  is  real  then  the  solution  of 


X 

(2.13)  mv  +^{v  f  ^[^jf  =  f(^) 


-1  <  4  <  1 


is 


)'''    '\ 


10 


2.1^)     iii)  =  iJl 


2      r  u- 


t'^)f(t)dt  ^2^,  r(|)  +  x2 


tTT 


+  ^l2A(k^+4k^)-k2i 


It  is  easy  to  verify  that  (2.12)  satisfies  (2,1)  when  S(z) 
satisfies  condition  (A).   If  we  write 


(0   =   2?i    f^o^^)  ■'^TTi   f(a] 


^[s^(a-2^i  f(c)] 


we  have 


r<b(T)dT  1         r^[S^(T)+27ri   f(T)]2_  [S^(T)-27ri   f(T)]2idT 

J  -T:rc=  271  J  T^ 

L  L 

r  (i)(T)dT  _   1    r/v  s(T)dT 

L  C 

Now  if  we  detach  C   from  L  and  expand  C   into    P   we  have 


/4^  =  |[[So(a+27rl   f(a]^+[s^(C)  -27ri   f(a] 


S(T)dT 


2?i  r  -TT-r 


and,    since  S(t)  — >A  as  t  — >■  oo 


y^iilJdT  ^   l|"|^S^(^)^    27ri   f(a]'^"+[So(^J-2^i   ^(^^F]-'^      • 


L 

Therefore 


11 


L 

-  2?i  [^''^^^  -S'(0 


The  condition  that  S(z)  must  be  analytic  in  the  complex  plane 
slit  along  L[-l,l],  [the  condition  (A)],  does  not  in  general 
determine  the  constants  k  and  k,  numerically.  This  can  be  seen 
by  considering  the  equation 

1 

(2.15)  Acl)(^)+Ml)  f    ^11^=   0  . 

J  L/  -  S 

-1 


The  function  S(z)  for  this  equation  is 


S{2)  = 


o   k  -  2Ak,  -  2Ak  z 
2  ,   o     1     o 


■^77 


'x^z^ 

-SAk^z 

< 

-  2Ak^  - 

-A^- 

z^- 

1 

and  this  certainly  satisfies  condition  (A)  if  vje  take  k  =  k,  =  0. 
For  these  values  we  have 


Ui)   =  |i  (A^)^  =  0 


namely  the  trivial  solution  of  (2.15).   On  the  other  hand,  S(z; 
will  satisfy  (A)  if 


12 


2  2  2  2 

A  z  -  2Xk  z  +  k^  -2Ak,  -  A^ 
o    o     1 


possesses  a  double  zero,  that  is,  if  vje  take  2k,  +  A  =  0.   With  this 


choice,  k  is  arbitrary,  and 


Az  -  k 
S{z)  =  '- 


from  which 

,  ft  A^ -  k^   ,   a4  -  k^ 

yr-i       Ji-r 

and  this  is  the  general  solution  of 

1 


x+l  i(t)dt^  0 


Thus  we  can  see  that  in  order  to  fix  k  and  k,  uniquely  we  may  have 
to  impose  side  conditions  in  addition  to  condition  (A). 

Let  us  examine  the  possibility  of  having  a  solution  ^   whose 
endpolnt  behavior  is  such  that 

1 

(B)  £^^^  (a-T)2  (()(t)  =  0  . 

With  this  behavior  the  endpoint  behavior  of  f(0  must  be  such  that 

(B^)  £^^^^  (a-af(a  =  0  . 

Under  condition  (B)  the  order  of  integration  can  be  changed  in 

|^(^)^r^(T)_dTdc 

L      L 
and  hence  by  integrating  (2.1)  we  find 


13 


L  L  L  L 

^  J ^{(Ddi:  +  f  miT)  -f(T)]dT  =y  fcodc 


(2.16) 


L 


Also,   after  multiplying   (2.1)   by  C  and  integrating: 

A/a(Odc  -/MT)/^^-^pJ^^^)^^^-i=/cf(adc 

L  L  L  L 

r  T  2 

Ajr^,|)(adC-      J(|)(T)dT         +y  t[A(|)(t)-  f(T)3dT   =  f  ^T{^)d^ 
L  [  L  J  L  L 

/ 


(2.17) 


2Ak^  -2   /    Tf(T)dT  =  k^    . 


The  substitution  of  these  values  for  k  and  k,  in  (2.9)  gives 

1 


(2.18) 


S(z)  = 


TfMdi 

J         T  -  Z 


and  a  solution  satisfying  (B)  will  exist  if  f(0  satisfies  (B^)  and 
(2.18)  satisfies  (A).  When  this  is  the  case  the  solution  is 


^'o  =  ^r2/£i^ 


27ri  f(C)  +A' 


1 


rf(T)dn 


2    I   ±-1Jl4^  _  2tt1   f(C)  +A^ 


14 


Instead  of  imposing  the  condition  (B)  at  each  endpoint  of  L 
we  may  wish  to  impose  it  at  only  one  endpoint,  say  the  condition 


(C) 


£^^1  (1-T)2  c1)(t)  =  0  . 


The  function  f(C)  would  then  have  to  be  such  that 
(C^)  £^.^1  (l-C)f(C)  =  0  . 

Under  condition  (C)  we  find  from  (2.1)  that 

>/(i+a<l>(adc  -/i(T)/i^^i±i+^M^M^  =  /(i+af(ad^ 

L  L      L  L 

1 

T^y  (i+a<l)(c)dc-   /<i)(T)dT  ^4-/" (i+T)[At|)(T)-f(T)]dT  =y' (l+af(Ud^ 

L  [  L       J    L  L 

(2.19)         2A(kQ+k^)-k^  =  2  /  (l+T)f(T)dT  . 

L 

The  substitution  of  this  value  of  k,  in  (2.9)  gives 

2Ak 


(2.20) 


S(z) 


2_   r  (l+T)f(T)dT   '-""o  ^  ^2 


1+2  J 
L 


1  +2 


and  a  solution  satisfying  (C)  will  exist  if  f(C)  satisfies  (C, )  and 
(2.20)  satisfies  (A).   When  this  is  the  case  the  solution  is 


'(^)  =  2FI 


r(i+T)f(T 

J   T^ 


)dT 


2Ak^    Q- 

+  2Tri  f(C)  -  YTT   +  ^ 


1 


2     r  (l+T)f(T)dT     ^  .  ^if..  2^^0 

+r  J  rr? 27ri  f(c)  -  YTT 


i+r 


+A' 


15 


3.  Equation  II 

In  this  section,  we  show  that 


(3.1) 


7r^i^(a  - 


J         T- 


=  f(a 


can  be  converted  into  a  particular  case  of  (2.1).   We  saw  in 
Section  2  that  if  (t5(T)  belongs  to  the  class  ^,  then  we  have  the 
identity 


(3.2)      2     r  ii^^Hui   r 


cl)(T)dTdC 


=  -A^(CD)  + 


/ 


tl)(T)dT 
T  -CD 


1  ^ 
1-CD^ 


/  (|)(T)dT 
L 


Hence  if  ^{r)    satisfies  (3.1)  then  it  must  satisfy 

(3.,)  ^  r  ii^^mn  rkis}^  =  .ri.)  .S 

or 

(5.4)  2/liz|!iiUi    lilildTdC  __   .,i.,2,f,,,    ^  ^2 


where 


^O    "  J     't^^'f^dT 


Equation  (3-^)  is  a  singular  integral  equation  of  the  first 
kind  which  can  be  solved  for 


2(l-r 


)Ma/i^ 


)dT 


16 


The   solution  is 

r 


(3.5)        2{l.^^)H0f^ 


L 

^                   nJl-T^     [k^    -(l-T2)f(T)]dT 
J     ^^ 


where 

L  L 


L 
.2 


L  L  L  L 

and  therefore 


TTk  =   f  Hr)  f  (C+T)<|>(adCdT 


L  L 


(3.6)  TTk    =    2k^k^ 


From  (3-5)  and  (3-6)  we  find  that  if  <|)(t)  satisfies  (3-1)  it  must 
satisfy 

(3.7)  MC)/i^=  F(C) 


where 


(3.8)   F(a 


1  Hl-T'^) 

2.2(1-^2)5/^  j     f:rr 


(l-T2)5/2^(T)dT 

"^rrx 

^^2      k  k 


'  2.(l.C^-)5/2  %(i.^^)3/2   • 
Conversely,  if  Mr)  satisfies  (3.7)  it  must  satisfy  (3-1)' 


17 


Equation  (3.7)  is  equation  (2.1)  with  X  =  0.   We  can  apply 
the  theory  of  Section  2.   It  follows  from  the  results  of  Section  2 
that  (3.1)  possesses  a  solution  provided 


S^(z) 


d-z") 


2  -1 

r  (l-T^)F(T)dT      ^O 

/  — ?ri r~2 


1-z 


satisfies  condition  (A),  i.e.,  S^(z)  is  analytic  for  z  not  on  L. 
Since 


2     r(i-'r^)F(T)dT  _     1        r       1         r(i-r)^/^f(adCdT 

r?  J        ^-^        ".^(i-z2)   J    ,^  ^,r— 2  i         ^--^ 


7r{ 


L   (t-z)/i-t-  L 

;^^ r  _rdT_^_fVi_   r      dT 

1--')  L  (T-z)fr7  -(1-')  L  (x-z)ji7;^ 


7r"(l-z") 
.2 


_i ^   r(i:£)£!f 


(C)dC 


ik^(2ki  +  k^z) 


d-z")       d-z") 


237F 


we  see  that  if  (3.1)  possesses  a  solution  then  the  constants  k^  and 
k^  must  be  such  that 


S^(z)  = 


1         ru: 


.c2)5/^f(C)dC  ^i^o(2ki  +  V^' 


X-z 


(l-z^)^/^ 


is  analytic  for  z  in  the  exterior  of  L.   If  this  condition  can  be 
satisfied  by  a  proper  choice  of  k^  and  k^,  the  solution  of  (3-1)  ii 
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(3.9) 


t)  = 


^71 


-     f(T)     + 


7r"(l-T") 


^^737^  / 


ik^(2k^+k^T) 


(1-t") 


3/2 


1 

27ri 


f(T)     - 


TT^d-T") 


27172 


/ 


(l-^^)^/^fU)dC 


C   -T 


lk^(2k^+k^T) 


For  the  case  in  which  L[-l,l]  coincides  with  the  real  axis 
we  find  that  if  f(C)  is  real  then  the  real  solution  of 


(3.10)     iT^^^ii) 


is 


(3.11)   <{>(?)  =  ^S 


r  (i)(t)dt 
J     t-g 


=  f(e)  , 


-1  <  e  <  1 


f(e)  + 


TTd-e") 


^737^ 


/ 


'  (l-t^)^^^f(t)dt 


TTT 


ik^2k^k^ 


provided  condition  (A)  is  satisfied. 

4.  Equation  III 

The  equation 


(4.1) 


7r¥(C)  + 


/ 


)(T)dT 


fca 


f,  .  .-Jj 
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can  be  written 


(^.2) 


'tti  Ma  +/iil:^].r-7ri  ^iO  -^f^^ 


no  . 


We  will  solve  this  equation  by  using  Carleman's  method.   Let  us 
introduce  the  function 


(4o) 


P(z) 


r<i)(T)dT 

J         T  -  Z 


a  function  of  the  complex  variable  z  =  x  + iy,  which  is  analytic  for 
z  not  on  L,  a  simple  smooth  arc  directed  from  z  =  -1  to  z  =  1.   If 
we  let  z  approach  a  point  ^  on  L  from  the  positive  side  of  L,   and 
then  the  negative  side  of  L;  the  limit  values  of  F(z)  are 


(4.4) 


(4.5) 


fUO   =  ^i  ct)(a  +/^^ 
L 

F-(c)  =  -TTi  Ma  +/4^ 


Equation  (4.2)  says  that  these  limit  values  must  satisfy  the 
barrier  equation 


(4.6) 


p+(a  •  F-(a  =  f(a 


where  C  is  not  an  endpoint  of  L.   The  function  P(z)  which  satisfies 
this  equation  and  has  the  properties  implied  by  the  representation 
(4.3)  provides  the  solution  of  the  integral  equation  because  the 
subtraction  of  the  Plemelj  formulas  (4.4)  and  (4.5)  gives 


(4.7) 


Ma  =  4t  tp^'K)  -F"(a]  . 


2TTi 


to?    Hi 
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If  we  write 

G(z) 


(4.8)  F(z)  =  ^ 


then  (4.6)  implies  that  G(z)  must  satisfy 
(4.9)  G"^(C)  +G"(C)  =  A  f(C)  +27ri  n 

where  n  is  an  integer.   As  can  be  seen  by  using  the  Plemelj 
formulas,  the  function 


(4.10)       G3_(z)  -^-^f    f    —^^i^Ihl^I-  +  ^1   n 


(jl-r^Yir-z) 


is  a  particular  solution  of  (4.9).  By  J  1-z^   we  mean  here,  and  in 
what  follows,  the  branch  of  the  two  valued  function  (1-z^)  '^'^^  which 
is  analytic  for  any  finite  z  not  on  L,  and  such  that 

Z  EH-      1 

Z->-00    z  ' 

The  function  G, (z)  is  analytic  for  z  not  on  L  and  since 

^z->oo^l^^^  "  const.  , 

^^  ^^^^  G^(z) 

Z  — =  0  . 

z-!^oo      I ^ 

z  +  LJ  1-z 
The  general  solution  of  (4.9)  is 

G(z)  =  G^(z)  +  /l-z^  p(z) 
where  p(z)  must  satisfy 


}-^ 
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(4.11)  P'^IO-P'IC)  =  0  . 

The  function  p(z)  must  be  taken  so  that  the  properties  of 


e^^^V  2  +ijl-z^  match  those  of 


This  function  is  analytic  for  z  not  on  L  and  it  vanishes  like  c  /z 

as  z  — >  00,  provided  c  =  /  (j)(T)dT  ^   0.   Furthermore,  in  accordance 

L 
with  the  assumptions  about  h^Kx)   admitted  in  the  introduction,  the 

behavior  of  F(z)  in  the  neighborhood  of  an  endpoint  a  of  L  is  such 

that  £  _^^(a-z)P(z)  =  0;  and  the  limit  values  P"''(C),  F"(0  must 

satisfy  a  uniform  Holder  condition.   These  properties  and  the 

condition  (4.11)  imply  that  p(z)  must  be  analytic  everywhere  and  it 

must  vanish  at  infinity,  i.e.,  p(z)  =  0. 

We  have  now  found  that 

,..is,  .<.,  = .  (..i/i:7>  exp  {lii!  /    ±y'^- 

L  Wi-t'^J(t-z) 
This  gives 


(4.15)  F^(c)  =  ±  (^-iiT^)yf(r7  .exp|%f.  \   .^iilzMi_| 

L  yi-T^(T-a  J 

(4.14)     F-(0   =  ±   k  .i/IIFj/fTTT    .  exp  \-   ^    f    AlLlMl.) 

^  L   ji-T2(T-aj 

The  solution  of  (4.1)  is  obtained  by  subtracting  (4.14)  from  (4.13). 
It  is 


,x^«  .= 


:.  jtv=,i::*;:nj<j"l  axfiT 


•fi   u;m  %t^v^ 


:vi:U     ^v.!ijjl    X?' 


^^ft 


■  m\. 
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(if. 15)   ^(O 


[F+(a  -p"(0] 


C  sin 


TT     \ 


2-rr  J       I o 


+  jl-C  cos 


L  /i-t'^(t-0 

^^^  f(T)dT 


/i-C 


2   r   ^..  ^/^^^J  I 
2Tr"  J 


L  yi-?'(T-ajJ 

^sis  we  have 
that  F(z)  =  /  ^^  _    vanishes  like  c^/z,  o^  ^  0,    as  z  — >  oo . 


It  will  be  noticed  that  in  the  above  analysis  we  have  assumed 

If 


this  is  not  the  case,  that  is,  if  for  example 

(4.16) 

while 

(4.17) 


r  (l)(T)dT  =  0 
L 


:^=-/xi( 


T)dT  =  0 


then  F(z)  vanishes  like  c^/z  ,  c^  ^  0,   as  z  — >  oo  and  the  solution 
of  (4.6)  as  we  have  given  it  has  to  be  adjusted.  However,  it  is 
clear  from  the  above  analysis  that  if  (4.l6)  and  (4.17)  prevail  then 
the  adjustment  is  easily  made  by  taking  F(z)  =  e   V(z+i^l-z  ) 
instead  of  (4.8),  which  leads  to 


4.18)   F(z)  =  ±  [z-i/ll?]^  expj  l^Ff-   / 


■^-^  f(T)dT 


(i^^y 


(T-Z) 
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5.  Some  Non-linear  Problems  in  Potential  Theory 

The  equations  vie  have  analyzed  can  be  identified  v.'ith  certain 
problems  in  potential  theory.  For  example,  consider  the  problem  of 
finding  ^(x,y)  such  that 

^xx^^'^^  ^^yy^^'^^  =  0  ,         y  <  0  , 
^y(x,0)  =  0  ,  |x|  >  1 

and  subject  to  an  additional  condition  on  y  =  0,  [x]  <  1  which  is 
specified  below.   Let  us  suppose  that  ij/   {x,0) ,    |x|  <  1,  satisfies 
the  conditions  imposed  on  (j)(t)  in  Section  1,  and  that  each  of 
f   (x,y)  and  ^  (x,y)  vanishes  as  z  =  x  + iy  — >  oo .   The  harmonic 
function  ^  (x,y)  can  then  be  written 

^  -'_^   (t-x)  +y 

1 

=  -h^  /^a(t-x)2+y2]^^(t,o)dt 

-1 

1 

^'   (x,y)  =  -  J-  f 


from  which 


•*■  (t-x)^  (t,0)dt 
•^   {t-x)2  +  y2 


and  by  integration 


1 

r 

_^   (t-x)''  +  y' 


,   r  (t-x)^  (t,0)dt 
^^(x,y)  =  i  /  ^ ^— 


r^    ^  (t,0)dt 

^x(->°)  =  I  j    \-X 
-1 


Now  if  we  impose  the  additional  condition 


2k 


Ayx,0)+^^(x,0)^y(x,0)   =Ii2L) 


x|    <   1 


then  (l)(t)   =  ^   [t,0),    \t\    <   1,   must   satisfy 


7.cl)(x)+(l)(x)    f    MtJ|t  =   f(x) 


-1 


If  instead  of  I  we  impose 


|x|  <  1 


II 


^l{x,o)-r^u,o)  =  lifl 


w\  <  1 


then  f(t)  =  ^  (t,0),  must  satisfy 


r  1 


7r^c|)2(x) 


r  (|)(t)dt 
J    t-  X 


f(x) 


1x1  <  1  . 


Finally,  if  the  additional  condition  is 


III 


^J(x,0)  +  ^^(x,0)  = 


1x1  <  1 


then  (|)(t)  =  ^  (t,0),  must  satisfy 


Tr^cp^{x)  + 


1       12 

r  (|)(t)dt 

J    t  -X 


-1 


f(x) 


1x1  <  1 


We  proceed  to  show  that  problem  III  can  be  solved  for  domains 
more  general  than  the  half  plane.   For  the  half  plane  problem  we 
can  write 

^(x,y)  =  (I  Jiz) 
and  then 
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9\z)   =  (£9'(2)  +1  ciJ^^'(z)   =  ^^(x,y)  -i^y(x,y) 


1  /o  r  (!)(t)dt  .  ±  'i  r  (i)(t)dt 


which  is  the  same  as 

(5.1)  .T/,.,=i/i(SJ|i    . 

-1 

The  integral  on  the  right  hand  side  of  (5.1)  has  been  determined  in 
connection  with  the  problem  of  solving  the  integral  equation  III, 
in  Section  h,    and  we  have 

(5.2)  ?^'(z)  =   ±   (^-^fl^)^exp[4#   /  -^^(t)at       1 

^^      -1    (il-t^J'^(t-2)  -^ 

where  y   is  1,  or  2.   The  function  f{:i,y)   can  be  obtained  from  (5.2) 
by  an  integration. 

Let  z  =  m(  C )  =  ni(  s  +iTi)  be  a  function  which  maps  the  domain 
D,  with  boundary  C  in  the  C-plane,  conformally  into  the  lower  half 
of  the  z-plane.   Let  the  image  of  C,  a  part  of  C,  be  the  segment 
y  =  0,  |x|  <  1;  and  let  the  image  of  Cg,  the  remaining  part  of  C, 
be  y  =  0,  |x|  =*  1.   Under  this  mapping,  M,(x,y)  is  transformed  into 

^(x,y)  =  tr(e,ri)  =  ^2  >  [m(a]  , 

a  function  harmonic  in  D. 

Let  s  be  the  arc  length  measured  along  C  from  say  the  initial 
point  of  C,  .   If  r=  4(s)  +iri(s)  is  a  point  on  C,    then  the  normal 
derivative  of  ^(s,!!)  at  the  boundary  is 
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(5.3)  -^3^r— =  lT^^^t-(^)5 

dm(.T)        dx 
nen  - 

Hence  for  cT  on  C 


If   cr  is  on  Cg,    then  ^^j^  =  if  ^^  ^®^^  ^"^   "J^'Lmlcr)]    is  real. 


2 

111   (4,Ti)   =   0 


n 


The  tangential  derivative  of  ^{^,r\)    at  the  boundary  is 


From  (5.5)  and  (5.4)  we  see  that  for  (Ton  C, 
— n       ~"s        I  as   I  I  1 


= 

dm{o') 
ds 

^[i'il-T[m( 

[■r 

= 

dm{cr) 
ds 

• 

Therefore 

if 

we 

take 

h( 

s)   =     ■ 

am(<r)  ^  f[ni(cr 
,2 

as  prescribed,  then 

i(e,ri)  =£  '?'[m(C)] 

is  harmonic  in  D,  with  normal  derivative  equal  to  zero  along  C^; 
and  with 
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ili^(e,n)  +11^(1, n)  =  h(s) 


prescribed  along  C,.   Thus  the  solution  of  problem  III  for  the 
domain  D  can  be  found  by  substituting  z  =  m( C ) ,  t  =  m(y)  and 


f(t)  =  f[m(cr)^  -  ^^h(^) 


dm(-T) 


in  (5.2).   The  substitution  gives 

(5.5)   d2u(oi.±-aci[,(a-ili:?[7)]' 

H    ds 


exp  J  -  il-LUI   r 

®^P  277  i  J 


'2fy^    r     ^^ll^^^^^Tr^l    •h{s)i-  ^^^g^  ds 


C^         Jl-m''^(cr)[m(o')  -  m(  O  ] 


and  then  1|J(4,ti)  can  be  obtained  by  integration. 

The  solution  of  problem  III  for  particular  domains  and  special 
boundary  values  of  the  gradient  is  not  new.   For  example,  it  is  vjell 
known  that  the  method  of  conformal  mapping  can  be  used  to  find  a 
potential  function  tC^^il)  such  that 

ii  riTi 


{^,r\)    =   0  T]    =    -1,1         ;    -oo    <    4    <   0      ; 


i^l(i,r])  +i/^{'i,r])   =   const.  n  =   -1,1        j        0  <   ^   <  oo    . 

This  problem,  or  its  equivalent,  occurs  in  the  hydrodynamical 
theory  of  Jets.   However,  the  author  has  not  seen  the  formula  (5.5) 
in  the  literature  and  it  is  probable  that  there  are  several  problems 
in  field  or  flow  theory  for  which  it  would  be  advantageous  to  use 
(5.5)  directly. 
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In  order  to  see  a  simple  application  of  (5.5)  consider  the 
problem  of  finding  jl(^,r))  such  that 

i,   (^,Tl)  +jj   (e,Tl)  =  0  ,       -00  <  I  <  00   ;    0  <  T)  <  1 

(5.6)  ilJ  (1,0)  =  0  ,         -00  <  ^  <  00 
t^Ai,l)  +fu,l)   =  h(|)  ,     -00  <  e  <  00 

The  function 

(5.7)  z  =  x  +  iy  =  coth  5  (l  +  iri)  =  coth  g  C  =  m(  C ) 

maps  the  infinite  strip  into  the  lower  half  of  the  3-plane  and  the 
segment  4'-2  =  0;  ji-^^lz)!  <lls  the  image  of  C,  :   C  =  l  +  i; 
-00  <  4  "  CO  .   From  (5-7)  we  find 


m'(a  =  -  ?  csch^  ^  C 


Jl-m^(0  =  v  1  -  coth^  ^  =  -i  csch  ^ 

m(o-)  =  m(A+i)  =  tanh  ^   A 
Jl-m^(r)  =  Jl  -  tanh^  ^  =  sech  ^  ^ 


dni(cr)  ^  _d_  ,     Z  A )  lil  ^  =  .1 

ds     dA  ^^^"^  2  ^^  ds   '  ds  . 


I  sech^  5  A   , 
cosh  ^  A  .  sinh  J  /; 


1     _    ""'^^^  2 
mUT)  -m(C) 


cosh  I  (A-  C) 


The  substitution  of  these  quantities  in  (5.5)  gives 
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(5.8)   f^  ?[m(s)]  =  ± 


2  smh  g  ^ 


cosh  K  C  -  1 


exp 


[i^ 


slnh  ^  r 


4  TT 


^  ^   i^-[4.cosh^  5  A.h(A)] 


dA 


cosh  ~   (A-O 


which  if  we  use 

,  1^  Jly.{h   cosh^  5  A]dA 

(5.9)  ^y 


-CD 


cosh  \   (A-C) 
=  ^2  +  2  £v»[l+cosh  I  C]  +27riTi 


becomes 


(5.10)     ^  3^[m(a] 


2   sinh 


•    exp     ^ 


.     ^2    TT 


'cosh  5  ^-1*''' 


/ 


smh  -^  u, 


h(A)dA 


.[cosh  5  C+l] 


-co    ^°^^   2    (^-^) 


If  we   choose  7  =   2   so  as   to  avoid  a   singularity  at   ^  =   0  we   finally 
have 


(5.11) 


d? 


■j^[m(C)] 


exp 


CO 


h(A)dA 


L"^'         cosh  ^  (A-O 


which  determines  the  derivatives  of  ^(^,t)).   That  is,  if 

t(e,Ti)  =  ^M(C) 
is  required  to  satisfy  (5.6)  then 


M'(0  =  ±  exp 


I-/ 


h(A)dA   ) 
cosh  ^  (A-O  i" 
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6.  Some  Non-linear  Systems 

In  a  recent  paper,  A.  Beurllng  [5]  proved  that  If  D  is  a 
region  in  Euclidean  n-space  then  the  non-linear  system 

(|)(x)  J  k{t,x)f{t)dt  =   f(x) 
(6.1)  D 

^(x)  y  k(x,t)<i)(t)dt  =  g(x) 
D 

possesses  a  unique  solution  when  the  kernel  k(t,x)  and  the 
prescribed  functions  f(x),  g(x)  are  suitably  restricted.   Since 
Beurling ' s  analysis  does  not  cover  the  following  singular  one- 
dimensional  case 


(6.2) 


^U)f  4ll|i=f(x) 


:x)/Mi) 


^(x)  /  :iii^  =  g(x) 


and  since  this  system  is  related  to  certain  problems  in  potential 
theory;  it  is  of  some  interest  to  see  that  (6.2)  can  be  reduced  to 
a  linear  system  by  using  the  method  of  Section  2. 

VJe  suppose  that  each  of  the  unknown  functions  belongs  to 
the  class  ^-j,,    defined  in  the  introduction,  and  we  suppose  that 
each  of  the  prescribed  functions  f(x),  g(x)  belongs  to  the  class 
0^2'      Then,  by  using  the  Hardy-Poincare-Bertrand  formula  as  shown 
in  Section  2  it  is  easy  to  verify  that 
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(6.3) 


where 


I'V 


f   (l-x2) 


irrr 


-1 


^(x)  /  4i)|t^^(^)|  iu)|t 


dx 


a_b. 


1-r 


=   -TT^^dWU   +/    4^At^-/    4^ 


^0=/    ^ 


(t)dt 


X 

^0=1     ^(t)dt 
'-1 


If  (6.2)  possesses  a  solution,  then  (6.2)  and  (6.3)  imply 
1 


(6.4)  -1^  [  (l-x-)[f(_x).g(x) 


1-1 


-1 


]dx  _  "o'^o 

1 


.,2^(i)^(e)  -./  4^-/  4_^ 


t)dt 


-1 


or  after  multiplying  by  ^{^)f{i) 


(6.5)  ^V(U^2(^)  .^iMl^ 

i-r 


r  (i-x^)[f(x)+g(x)]dx       ^ 

7  X  -  ^  o  o 


-1 


.  ^(^)/  4t)dt.^(^)J  iU)dt^  0 , 


-1         -1 

Using  (6.2)  again,  equation  (6.5)  becomes 

(6.6)        ir^\^^(i)^^{i)   +iiiMil.p(e)_f(^)g(4)  =  0 
where 


1-^ 


F(|) 


-I^ 


x'^)[f(x)  -hg(x)]dx 

— 5rri 


-  a  b^  . 

o  o 
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Equation  {6.6)   fixes  the  product  ^{^)f{i).      We  find 

(6.7)    i(^we)  =  m)  -  -  miil^'ii^-^f^i^Mi'on-eji . 

With  this  product  knovm  the  system  (6.2)  can  be  written 

1 
r    ^(t)dt  _  f(x)  ,   ,  . 

(6.8) 

1 


/  -¥-¥  =  fin  •  *(>=>  • 


The  last  equations  can  now  be  solved  by  Carleman's  method  [^]. 

Notice  that  if  each  of  f(x)  and  g(x)  is  integrable  and  we 
require  each  of  the  unknown  functions  to  satisfy  the  condition  (B) 
of  Section  2  then 


1  11  11 

dxdt 
t 
1         -1      -1  -1      -1 


^1  -1 

1 

=  -J    g(t)dt 

-1 

or 

1 
(6.9)  /  [f(x)  +g(x)]dx  =  0  . 


Also, 

1  11 

J    xf(x)dx  =  ~f    fit)  f    iiE^tttl  <|)(x)dxdt 


-1    -1 

-/  ^(t)dt./  <|>(x)dx  -/  tV/(t)/  il|)^ 

-1      -1      '-1     -1 
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which  determines  a^b^,  that  is. 


1 

h^  =   f    x[f(x)  +  g(x)]dx  . 


0  0    J 

-1 

Hence  for  solutions  which  satisfy  (B),  F(^)/(l-|")  becomes 


F(i)  ^  r  [f(x)+g(x)]dx  __g     r  iru)^gu)] 
i-r    ^1       ""'^  i-r  ^_3^ 


dx 


and  since  the  second  term  on  the  right  is  zero  by  (6.9)  we  have 

1 


Fi'c.)  ^  r  [f(x)+g(x 


)]dx 


.1  <  X  <  1 


We  remark  in  passing  that  the  system 
f^   VT(t)dt 

-1 

p    <|)-,(t)dt 
-1 
can  be  reduced  to  (6.2)  by  using  the  substitutions 

<j)^(t)  =  M-t)   ,    f3^(x)  =  f(-x)   , 
^3^(t)  =  ^(t)    ,    g^(x)  =  -g(x)   . 

The  systems  noted  above  are  not,  of  course,  the  only  non- 
linear systems  that  can  be  solved  by  using  the  technique  of 
Section  2.   The  Hardy-Poincare-Bertrand  formula  can  be  used  to 
linearize  a  number  of  other  systems  in  v/hich  the  non-linearity  is 
due  to  the  presence  of  a  product  in  which  an  unknown  function  is 


q   O/j.-'J"    -^v    f'lib 


3^ 


multiplied  by  an  unknovm  Hilbert  transform  or  an  unknown  Stieltjes 
transform.   As  two  additional  examples  we  cite  the  systems 

1 


and 


-1 

1 

f{x)f    %^-|^=  g(x)+A2^(x) 
-1 

0  0 


The  last  system  plays  an  important  role  in  the  theory  of  radiative 
transfer. 
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